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The energy per particle of many body wavefunctions that mix Laughlin liquid with crystalline 
correlations for periodic samples in the Haldane-Rezayi configuration is numerically evaluated for 
periodic samples. The Monte Carlo algorithm is employed and the wave functions are constructed in 
such a way that have the same zeroes as the periodic Laughlin states. Results with up to 16 particles 
show that these trial wavefunctions have lower energy than the periodic Laughlin states for finite 
samples even at v = | . Preliminary results for 36 particles suggest that this tendency could reach the 
thermodynamic limit. These results get relevance in view of the very recent experimental measures 
that indicate the presence of periodic structures in the 2DEG for extremely small temperatures and 
clean samples, inclusive at main FQHE filling fractions v = |, |. 
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Although the Fractional Quantum Hall Effect is de- 
scribed in its relevant aspects in terms of Jastrow-like 
many-body wave functions, a link between this varia- 
tional approach and hamiltonian solutions is still desir- 
able. The states introduced by Laughlin cleanly incorpo- 
rate the tendency of particles to be as far away as possi- 
ble from each other. On the other hand the Hartree-Fock 
states discussed in the literature (See pj and the refer- 
ences therein) are approximate solutions of the Hamilto- 
nian and may have valuable information not taken care of 
by the variational states. With this possibility in mind 
we have constructed a new type of trial wave function 
that incorporates a part of each approach as described 
below. Precisely in these days this construction gets rel- 
evance thanks to very recent experiments that had de- 
tected the presence of surprising periodic structures in 
the 2DEG for highly clean samples subject to extremely 
small temperatures [2] . Assuming that the proposed here 
wavefunctions maintain their properties of showing lower 
energies than the so called periodic Laughlin states (in 
the Haldane-Rezayi periodic scheme Q) in the thermo- 
dynamic limit, they could have the opportunity of de- 
scribing the detected periodic structures. Qualitatively, 
the energy dependence on the sample size obtained here 
suggests that the energy is lower in the limit of infinite 



size samples. However, this important conclusion needs 
for more extensive calculations for its full confirmation, 
that are expected to be considered elsewhere. 

Analytic Hartree-Fock solutions in the lowest Landau 
level may be written down for filling fractions of the form 
v = 1/q. This unusual property was discovered after 
numerical results showed that there is a self-consistent 
charge density wave solution (CDW) to the Hartee-Fock 
(HF) equations that has a zero of order q— 1 in each CDW 
plaquette. This property was then used to diagonalizc 
the eigenvalue matrix. The resulting charge density cor- 
responds to a lattice of holes, with percolating ridges that 
surround the zeroes. The form of the associated single 
particle wavefunctions suggests a way of constructing our 
new state incorporating both the Laughlin correlation in 
part, as well as the HF crystalline correlations(See [H). 

The appropriate framework to achieve this on the plane 
using numerical procedures is to impose periodic bound- 
ary conditions on the single particle states, as suggested 
by Haldane and Rezayi many years ago [3| . This is done 
in Ref. The main idea is to exploit the fact that 

the mean field Slater determinant in this picture can be 
written as the product of another determinantal function 
containing the whole dependence on the quantum num- 
bers of the single particle HF states, times a factor whose 
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zeroes are spatially fixed and periodic with the periodic- 
ity of the density. Therefore, the position of those zeroes 
has no dependence whatever on the set of quantum num- 
bers of the filled mean field states, and has thus been 
factored out of the Slater determinant. Moreover, the 
number of zeroes of those kinematical factors as a func- 
tion of any of the identical particle coordinates is just 
(q — l)N e , where N e is the number of electrons. 

Let us begin by writing the explicit form of the counter- 
parts of the Laughlin wavefunction in the Haldane-Rczayi 
scheme for implementing periodic boundary conditions 

* i= exp(- £ lL){ M l(Z*-R*)\-T*)} q x 

i=l,2...JV e ° 

n {mI(z*-z*)\-t*)}\ 

i<j 
i,j=l,2...N e 

where z* = Xj —i yj , j = 1, 2, ....N e and Z* = X — iY are 
particles and center of mass complex coordinates, respec- 
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tively. L is the cell size, R* — niai+n 2 a 2 with ri\ and n 2 
integers, t* = exp{— 2iri/Q), and a\ = ^J2i:q/ sin(27r/6), 
a 2 = a\T. The functions -&\{u\t) are the odd elliptic 
theta functions and vanish as the first power in u as this 
variable goes to zero. We notice that as one particle 
approaches another, vanishes as a power q, so that 
altogether it includes q(N e — 1) zeroes of this kind, plus 
q generated by the center of mass factor. We can thus 
replace the factors with spatially fixed zeroes in the HF 
solution, by a proper Laughlin factor to obtain the same 
short range behavior (a zero of order q) when any two 
particles approach each another. However, the presence 
now of the detcrminantal function keeps the crystalline 
information associated with the optimization of the mean 
field problem. Therefore, the proposed states have an a 
priori chance of lowering the energy per particle of the 
Laughlin states. 

Then, in the Landau gauge A = — B(y, 0,0), as it is 
proper of the Haldane-Rezayi scheme, the state ansatz 
being proposed here has the explicit form 



i=l,2...JV e ° 
1 l^ 1 1-r (IT 1 l^ 1 

-0} n {M^t-m-r*)} , 

i,j = l,2...Ne 

(1) 



where Q = — ^ (q — 1),Z = Y^j=i z i an d r o = \Jhcj\eB\. 

The functions navc tnc f° rnl 1 1 

si = 5-11 x a 2 , s 2 = — x ai 

xL 0) (x) = exp(ix z* |^)x HMj;^ -R*- C k )| - r*). n = (0,0, 1), a r Sj = 2tt%, 

while the allowed values of n\ and n 2 are 



Here x = — k.ai/ai and the argument Ck depending on 
the quantum number k is given by 

C k = — (-k.a 2 - k.air ) H — . 

Ztt Z 

Our trial wave function "J is periodic, with a slanted pe- 
riodicity region of equal sides length L = Nai, and a 
slant angle of 2n/6. The number of particles in the re- 
gion is N e — N 2 . The momenta allowed by the periodic 
boundary conditions are 

, «i n 2 

k = 7/ Sl + 7/ S2 ' 

where the reciprocal lattice unit vectors and the normal 
vector are defined by 



ni + l,..,0,...,y-l}, 

n 2 e{--,-- + l,.,0,...,y-l}. 

The evaluations of the energy per particle of both 
states were done by employing the Monte-Carlo method 
for samples having a number of particles N e equal to 
4 and 16 for the case of the Laughlin states. As for 
the trial wavefunctions investigated here the calculations 
were done for 4, 16 and 36 particles. The results are illus- 
trated in the Tables I and II and in Fig. 1. As mentioned 
before, the parameter £ is the one defining the probability 
of admission of new configurations as usually is needed to 
do in the Monte-Carlo algorithm. The margins of errors 
reported correspond to the maximum deviation from the 
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TABLE I: The results of the energy per particle for the coun- 
terparts of the Laughlin states in the Haldane-Rezayi periodic 
boundary conditions. The state for JV e = 4 was evaluated two 
times for different values of the Monte Carlo new configura- 
tion admission coefficient £ in order to check the independence 
of its value. The calculations were done for N e — 4, 16 parti- 
cles 



N e = N 2 


€ 


e 


4 


0.3 


-0.374476 ± 0.0000697292 


4 


0.25 


-0.374245 ± 0.0000779082 


16 


0.3 


-0.392032 ± 0.0000880513 



TABLE II: The results for the energy per particle for the 
trial state proposed in this work. The same Haldane-Rezayi 
periodic boundary conditions were employed. In this case the 
states for N e = 4, 16 were evaluated two times each one for 
different values of £ to check the independence of the result 
on this constant. The evaluations were done for N e = 4, 16 
and 36 particles in this case 



N e = N 2 




e 


4 


0.2 


-0.414112 ± 0.0000625317 


4 


0.25 


-0.414191 ± 0.0000457599 


16 


0.25 


-0.410156 ±0.00006590 


16 


0.2 


-0.410222 ±0.00008659 


36 


0.25 


-0.410943 ± 0.000486141 



mean value of a set of the last 60 percent of the evaluated 
energies in the Monte-Carlo iterative process. 

The expectation value of the many-particle Hamilto- 
nian in the Laughlin state was evaluated using the Monte 
Carlo method for samples with N e = 4, 16. For our trial 
wavefunctions calculations were done for 4, 16 and 36 
particles. Results are shown in Tables I and II. The pa- 
rameter £ defines the probability of admission of new 
configurations in the Monte-Carlo algorithm. The errors 
reported in the third column correspond to the mean 
square of the fluctuations in the Monte Carlo output af- 
ter convergence was assured. 

Our results are plotted in Fig. 1. The error bars are not 
resolved at the scale of the plot. 

Note that a linear extrapolation to the thcrmodynam- 
ical limit N— > oo of the Laughlin state energy repro- 
duces former estimates of the energy of this state at 
v = 1/3 . It is also clear that the energy per parti- 
cle of each of the computed Laughlin values lies above 
that obtained for our trial wave function for any of the 
evaluated finite size samples. Moreover, the behavior of 
the energy of the latter for the largest number of parti- 
cles evaluated, suggests that the energy per particle in 
the thermodynamic limit is lower than the one associ- 
ated with the Laughlin state. In the contrary case, the 
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FIG. 1: Energy per particle as a function of the inverse square 
root of the number of particles N = y/N e . The picture shows 
that the introduction of correlations in the HF crystalline 
states made their energies lower than the ones shown by the 
versions of the Laughlin state in the Haldane-Rezayi periodic 
scheme. Note also that, if the behavior in the large N limit 
is confirmed by more extensive evaluations, the results will 
imply the existence of a ground state with a slightly lower 
energy than the Laughlin one for macroscopic samples. The 
curve joining the points of evaluated energies for the new trial 
state is a fitting of these three points to a quadratic polyno- 
mial in 1/ y/N e . The lower straight line with negative slope is 
simply a linear curve of l/y/N e minimizing the mean square 
deviations from the three measured energies. 

extrapolation curve in the variable 1/N which join the 
there three evaluated points of the energy N — 2,4 and 
6 should change its monotonic change of slope upon en- 
larging the value of N. However, being the number of 
Monte Carlo method iterations for the 36 particles state 
(N = 6) yet limited, this indication is not yet conclusive 
and further numerical evaluations will be done to give a 
better foundation to this conclusion. Its validity, clearly 
leads to the idea about that the recently detected peri- 
odic structures in extremely perfect 2DEG at very low 
temperatures, could be associated to the here proposed 
translation symmetry breaking states 

In ending we would like to underline that the present 
letter consider a particular HF state showing one electron 
per its periodicity unit cell. For fillings of the form 1/q 
the HF solution produces a gap for all integers q. The 
experiment suggests, however, that even and odd q values 
are qualitatively different states. It has been shown in 
the past that if only half electron is captured by the unit 
crystalline cell this distinction is properly borne out 0. 
Future work will extend to cover such states, and shall 
be reported elsewhere. 
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